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The existence, stability and other dynamical properties of a new type of multi-dimensional (2D or 
3D) solitons supported by a transverse low-dimensional (ID or 2D, respectively) periodic potential 
in the nonlinear Schrodinger equation with the self-defocusing cubic nonlinearity are studied. The 
equation describes propagation of light in a medium with normal group- velocity dispersion (GVD). 
Strictly speaking, solitons cannot exist in the model, as its spectrum does not support a true 
bandgap. Nevertheless, the variational approximation (VA) and numerical computations reveal 
stable solutions that seem as completely localized ones, an explanation to which is given. The 
solutions are of the gap-soliton type in the transverse direction(s), in which the periodic potential 
acts in combination with the diffraction and self-defocusing nonlinearity. Simultaneously, in the 
^ I longitudinal (temporal) direction these are ordinary solitons, supported by the balance of the normal 

O ■ GVD and defocusing nonlinearity. Stability of the solitons is predicted by the VA, and corroborated 

' by direct simulations. 

^ ■ I. INTRODUCTION 

Recently, a variety of two- and three-dimensional (2D and 3D) solitons have been investigated in models based on 
the nonlinear Schrodinger (NLS) or Gross-Pitaevskii (GP) equations with a spatially periodic potential and cubic 
P5 nonlinearity, see a review [ij . The physical models of this type emerge in the context of Bose-Einstein condensation 
(BEG) 1^ yilM B 0: where the periodic potential is created as an optical lattice (OL), i.e., interference pattern 
^ formed by coherent beams illuminating the condensate, and in nonlinear optics, where similar models apply to 
photonic crystals @. A different but allied setting is provided by a cylindrical OL ("Bessel lattice"), which can 
. also support stable 2D and 3D ^Td\ solitons. Additionally, models combining a periodic lattice potential and 
^ ' saturable nonlinearity give rise to 2D solitons, that were predicted in Ref. HjII and observed in several experiments 
in photorefractive media, including fundamental solitons 12] and vortices 13]. It is also relevant to mention that 
lO ' experimental observation of spatiotemporal self- focusin g of light in silica waveguide arrays, in the region of anomalous 

', group- velocity dispersion (GVD), was reported in Ref. . 
' I ' In models with the cubic nonlinearity, these solutions were investigated in a quasi-analytical form, which combines 
(y-^ ] the variational approximation (VA) to predict the shape of the solitons, and the Vakhitov-Kolokolov (VK) criterion 
1 |l6l ] to examine their stability. Final results were provided by numerical methods, relying upon direct simulations 
' of the underlying NLS/GP equations. A conclusion obtained by means of these methods is that, unlike their ID 
. counterparts, multi-dimensional solitons in periodic potentials can exist only in a limited domain of the ( A'^, e) plane, 

where N and e are the norm of the solution and strength of the OL potential, respectively. The most essential 
. . limitation on the existence domain of 2D solitons is that N cannot be too small (in a general form, a minimum value 
^ of the norm, as a necessary condition for the existence of 2D solitons supported by lattice potentials, was discussed 
k>( in Ref. ^3)- Unlike it, e may be arbitrarily small, as even at e = the 2D NLS equation has a commonly known 
s_j weakly unstable solution in the form of the Townes soliton, at a single value of the norm, AT = iVx fll (A^T ~ 11-7 for 
d ^ the NLS equation in the usual 2D form, iut + V^u -I- \u\'^u = 0). Small finite e gives rise to a narrow stability region, 

< iVx - iV < (AiV)„_ ~ e (1) 

for the 2D solitons ^ . Grossing the lower border of the existence domain |^ leads to disintegration of the localized 
state into linear Bloch waves (radiation) |0 . 

In the case of the attractive cubic nonlinearity (which corresponds to BEG where atomic collisions are characterized 
by a negative scattering length, while this is the case of the normal, self- focusing Kerr effect), 2D and 3D solitons 
can be stabilized not only by the potential lattice whose dimension is equal to that of the ambient space, but also 
by low- dimensional periodic potentials, whose dimension is smaller by one, i.e., 2D and 3D solitons can be stabilized 
by a quasi-lD 0, 0] or quasi-2D J_, _6, Tj OL, respectively [in the former case, the qualitative estimate ^ for the 
width of the stability region at small e is correct too]; however, 3D solitons cannot be stabilized by a quasi-lD lattice 
potential |^ Q [this is possible if the ID potential is applied in combination with the Feshbach-resonance management, 
i.e., periodic reversal of the sign of the nonlinearity coefficient 20], or in combination with dispersion management, 
i.e., periodically alternating sign of the local GVD coefficient |21|]. Solitons can exist in such settings because the 
attractive nonlinearity provides for stable self-localization of the wave function in the free direction (one in which the 
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low-dimensional potential does not act), essentially the same way as in the ID NLS equation, and, simultaneously, 
the lattice stabilizes the soliton in the other directions (in the 3D model with the quasi-lD OL potential, the self- 
localization in the transverse 2D subspace, where the potential does not act, is possible too, but the resulting soliton 
is unstable, the same way as the above-mentioned Townes soliton). An important aspect of settings based on the 
low-dimensional OL potentials is mobility of the solitons along the free direction, which opens the way to study 
collisions between solitons and related dynamical effects 0. 

In the case of defocusing nonlinearity, which corresponds to a positive scattering length in the BEG, or self-defocusing 
nonlinearity in optics (negative Kerr effect), the soliton cannot support itself in the free direction. Localization in 
that direction may be provided by an additional external confining potential; however, the resulting pulse is not a true 
multidimensional soliton, but rather a combination of a gap soliton (a weakly localized state created by the interplay 
of the repulsive nonlinearity and periodic potential 2]), which was recently created experimentally in a ID BEG 22]) 
in the direction(s) affected by the OL, and of a Thomas- Fermi state, directly confined by the external potential in 
the remaining direction 

Thus, no soliton can be supported by a low-dimensional lattice in the BEG model (GP equation) with self-repulsion 
(the latter corresponds to the most common situation in the experiment '2^1). On the other hand, a new possibility 
may be considered in terms of nonlinear optics. Indeed, one may combine three physically relevant ingredients, viz., (i) 
an effective periodic potential in the transverse direction(s), while the medium is uniform in the propagation direction, 
(ii) self-defocusing nonlinearity, and (iii) normal GVD. The latter is readily available, as most optical materials feature 
normal GVD, in compliance with its name. As concerns the negative cubic nonlinearity, it is possible in semiconductor 
waveguides, or may be engineered artificially, throu gh t he cascading mechanism, in a quadratically nonlinear medium 
with a proper longitudinal quasi-phase-matching |23 |. Also quite encouraging for the study of multidimensional 
solitons proposed in this work are recent observations of ID [23| and 2D solitons in optically induced waveguide 
arrays (photonic lattices) with self-defocusing nonlinearity. 

The setting outlined above can be realized in both 2D and 3D geometry, where the necessary transverse modulation 
of the refractive index is provided, respectively, by the transverse structure in a planar photonic-crystal waveguide, 
or in a photonic-crystal fiber. To the best of our knowledge, in either case the model is a novel one. A soliton 
in this medium, if it exists, will be of a mixed type: in the transverse direction(s), it is, essentially, a ID or 2D 
spatial gap soliton, supported by the combination of the effective periodic potential and self-defocusing nonlinearity, 
while in the longitudinal direction it is a temporal soliton of the ordinary type, which is easily sustained by the joint 
action of the self-defocusing nonlinearity and normal GVD. Thus, one may anticipate stable spatiotemporal solitons, 
alias "light bullets", in this model. Due to their mixed character, they may be called semi-gap solitons. The issue 
is of considerable interest in view of the lack of success in experiments aimed at the creation of "bullets" in more 
traditional nonlinear-optical settings . The only earlier proposed scheme for the stabilization of 2D spatiotemporal 
optical solitons in periodic structures, that we are aware of, assumed the use of a planar waveguide with constant 
self- focusing nonlinearity and longitudinal dispersion management l2^ . 

On the other hand, it is necessary to stress that, rigorously speaking, completely localized solutions cannot exist 
in the present model: its linear spectrum cannot give rise to any true bandgap, in which genuine solitons could be 
found (see below); instead, one may expect the existence of quasi-solitons, consisting of a well-localized "body" and 
small nonvanishing "tails" attached to it. Nevertheless, we will produce families of solutions which seem as stable 
perfectly localized objects. This is possible because the "tails" may readily turn out to be so tiny that they remain 
completely invisible in numerical results (possibly being smaller than the error of the numerical scheme), and, of 
course, they will be invisible in any real experiment. An explanation to this feature is provided by the fact that 
bandgaps, which "almost exist" in the system's spectrum, do not exist in the strict sense because they are covered by 
linear modes with very large wavenumbers. As shown in Ref. p^ . in this case the amplitude of the above-mentioned 
tails (which are composed of the linear modes with very large wavenumbers) is exponentially small. In fact, families 
of stable "practically existing" solitons in a second-harmonic-generating system with opposite signs of the GVD at the 
fundamental-frequency and second harmonics, where solitons cannot exist in the rigorous mathematical sense, were 
explicitly found in that system, in both multi- 26] and one- dimensional settings. Implicitly (without discussion 
of this issue), "practically existin g" s olitons (although, in this case, they were unstable against small perturbations) 
were also found in a recent work [28j, which was dealing with a 2D model of a planar nonlinear waveguide with the 
cubic nonlinearity, that features a Bragg grating in the longitudinal direction, and is uniform along the transverse 
coordinate. In the latter model, true solitons cannot exist, as the spectrum of the system does not support a full 
bandgap. 

The objective of the present paper is to explore 2D and 3D spatiotemporal solitons (which may be, strictly speaking, 
"quasi-solitons" , in the above sense, but feature completely localized shapes) and their stability in the proposed 
medium. In Section 2 we fix the mathematical form of the 2D version of the model, analyze its spectrum, and apply 
the VA to the study of solitons. In Section 3, direct numerical results demonstrating the existence of very robust 
2D solitons, and their delocalization when the lattice strength e becomes too small, are reported (the comparison 
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with the VA prediction shows that the VA provides for a crude approximation in the present model). In Section 4, 
we additionally consider the effect of variation of the nonlinearity coefficient along the propagation distance on 2D 
solitons. Finally, numerical results for 3D solitons are collected in Section 5, and the paper is concluded by Section 6. 



II. FORMULATION AND VARIATIONAL ANALYSIS OF THE TWO-DIMENSIONAL MODEL 

The model of a planar waveguide corresponding to the outline given above is based on the following variant of the 
2D NLS equation for the local amplitude u{z,x,t) of the electromagnetic field: 

iuz - + ^u^x + £ cos(2a;)u - = 0, (2) 

where z is the propagation distance, x is the transverse coordinate, and t is the reduced time, defined the same way as 
in fiber optics. The signs in front of the GVD {utt) and cubic terms correspond, as said above, to the normal GVD and 
self-defocusing nonlinearity, e is the amplitude of the transverse modulation of the refractive index (which is assumed 
sinusoidal, but the results will be nearly the same for more realistic forms of the modulation which correspond to the 
actual photonic-crystal structure), and the period of the modulation is scaled to be tt. Dynamical invariants of Eq. 
are the norm of the solution (in optics, it is the total energy), 

/ + 00 /' + 00 

dx / dt\u{z,x,t)\^, (3) 
-oo J —oc 

together with the longitudinal momentum and Hamiltonian, 

p+oo r + co 

P ^ i I dx I dt u*u, (4) 



-oo J —oo 

f+oo r+oo 

H = I dx dt 



OO 



i|u|''-ecos(2x)|up 



(5) 



To find the linear spectrum of the model, one looks for solutions to the linearized version of Eq. lO as 

u{z,t, x) — exp (ikz — icut) Fe{x), (6) 
where is a real propagation constant, uj is an arbitrary real eigenvalue, 

E = uj^-2k, (7) 

and Fe{x) is a solution of the Mathieu equation, 

i^" + 2ecos(2a;)i^ + £;i^ = 0, (8) 

corresponding to the eigenvalue E. It is commonly known that the Mathieu equation gives rise to bandgaps in its 
own spectrum, i.e., to forbidden intervals of the values of E, within which no regular quasi-periodic solutions of Eq. 
© can be found. However, since all large values oi E {E ^ e) belong to the allowed band, where such solutions exist, 
it is obvious that no value of k may fall in a forbidden bandgap. Indeed, using Eq. Q, one can construct any real 
value of the propagation constant as fc = (w^ — i?) /2, taking very large E and, accordingly, very large w. 

On the other hand, the same consideration suggests that, in some cases, the necessary values of E and u may be 
very large indeed. It was shown, in a general form, in Ref. |26j that short-period waves corresponding to such large 
parameters, which build up into a possible "tail" attached to the soliton's "body" (that makes it a quasi-soliton), 
will have an exponentially small amplitude, rendering the tail totally negligible (in particular, it may be completely 
invisible in numerical solutions). Therefore, it makes sense to look for "practically existing" solitons in the present 
model. 

We start searching for stationary solutions of the "mixed" type, which, as explained above, are expected to feature a 
gap-soliton (weakly localized) shape along x and strong ordinary localization in <, by adopting the following variational 
ansatz, 

u{z, X, t) = Ae'-^'' [x~^ sin(ax)] sech(at), (9) 
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where a and a are, respectively, the transverse and longitudinal inverse widths of the soliton, and A is its amplitude. 
Using the obvious Lagrangian representation of Eq. Q and well-known VA formalism , one can readily derive the 
following equations for the parameters of the ansatz, 



3e ^ 



1 + (iV/37r)^ 



(10) 



= K = £ r— — -, (11) 

Stt V67r 2 (277r2)i/3 (9^2 + 5Ar2)2/3 

where N is the norm defined by Eq. It follows from these equations that the condition dk/dN < 0, i.e., the 

necessary stability condition, according to the above - mentioned VK criterion p^ . always holds, as shown in Fig. ^ 
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FIG. 1: The propagation constant vs. norm for the 2D solitons, as predicted by the variational ap- 
proximation, Eq. Hll() . for different values of the strength of the periodic potential e. The nega- 
tive slope, dk/dN < 0, implies stability of the solitons according to the Vakhitov-Kolokolov criterion. 



The stability of the soliton families, predicted by the VK criterion as per Fig. ^ is, generally, corroborated by 
direct numerical simulations, see the next section (with a caveat that the VA predicts the shape of the solitons in 
only a qualitatively correct form, as explained below). However, it should also be mentioned that the applicability of 
the VK criterion to gap solitons in lattice models has never been proven, therefore one should apply this method with 
due care. In particular, it is known that the gap solitons may be unstable in some cases when they are expected to 
be VK-stable 30], which is not very surprising, as the VK criterion ignores complex stability eigenvalues, that may 
give rise to oscillatory instabilities. More perplexing is the fact that some gap-soliton families in a lattice model with 
the cubic-quintic nonlinearity, that are formally predicted to be VK-unstable, are in reality completely stable |3lj| . 

To complete the discussion of the VA, it is necessary to notice that ansatz @ is irrelevant if it predicts a ^ 1, as 
it would mean that the soliton is very broad in the x-direction, and it does not feel the underlying lattice structure, 
cos(2a;) in Eq. Wc therefore limit the applicability of the VA by a (roughly defined) condition, a > 1. According 
to Eq. (|l()|l . this implies that the potential must be strong enough. 



III. NUMERICAL RESULTS FOR TWO-DIMENSIONAL SOLITONS 

Direct simulations of Eq. jSJ (propagation in z) started with an initial localized waveform, which we took as ansatz 
with the parameters predicted by Eqs. Ulljl . or just a Gaussian with rather arbitrary parameters - for instance. 



u{x,t,0) = Aexp[--(x2 -t-t2)]. 



(13) 
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It was observed that the initial waveform undergoes intense evolution, shedding off some radiation waves that were 
absorbed at edges of the integration domain. The domain was large enough - in most cases, (— Stt, +87r) in both 
directions {x and t) - so that the solitons, which are typically well localized within a region (—5, +5), see Figs. 2, 3, 
7, and 9 below, are not affected by the absorbers. 



T 1 1 ' r 
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FIG. 2: Cross sections of a typical established 2D soliton of the semi-gap type. It has self-trapped 
from the initial configuration H13|) with A = I and a = 1. The solid and dashed lines rep- 
resent, respectively, the sections along the transverse directions x at t = 0, and temporal direction 
t at X = 0. The shapes of the cross sections are displayed on the logarithmic scale, to illus- 
trate the fundamental observation of the practical vanishing of the tails attached to the soliton's "body" . 

The evolution of the pulse ends with the establishment of a 2D stationary soliton, without any visible tails, as shown 
on the logarithmic scale in Fig. 2. It is relevant to mention that, as seen from Figs. 13 and 13 the characteristic sizes 
of the soliton in the x and t directions arc on the same order of magnitude, hence the values of the GVD coefficient 
and its effective lattice-diffraction counterpart are, roughly, equal. 

In fact, the solitons self-trap even from the initial configurations that are quite different from their final shape, 
which attests to strong robustness of the solitons. The soliton's stability was then additionally tested against small 
random perturbations, by simulating the evolution of the initial configuration U{x,t) = Ui){x,t)[l + (TUp{x,t)], where 
Uo{x, t) is the numerically found soliton, cr is a small amplitude of the perturbation, and Up{x, t) is a random function. 
An example of a stable localized state self-trapped from the initial Gaussian (|13|1 with A — 1 and a ~ 1, with the 
norm No = tt, is displayed in Fig. 13 

It should be said that direct comparison of the VA predictions with the numerical results shows only a qualitative 
agreement: for example, the VA predicts, with the same value of the norm, a soliton whose width in the temporal 
direction exceeds the actual width of the soliton in Fig. |2| by factor in excess of 2, and, accordingly, the amplitude 
of the numerically found soliton significantly exceeds that predicted by the VA. Therefore, the VA provides for only 
a crude approximation in this model; nevertheless, its qualitative predictions, such as the stability of the solitons 
predicted by the VK criterion, are correct. In this connection, it is relevant to mention that no good version of VA 
has been thus far proposed for gap solitons (this technical issue was considered, in some detail, in Ref. 29] ), and the 
solitons in the present model are still more complex objects. 

Further numerical analysis of the 2D model has demonstrated that, following the known pattern of the delocalization 
transition of 2D solitons in lattice potentials , the soliton solutions cease to exist when the strength of the periodic 
potential e or the norm of the localized state N fall below some critical values. In this case, the localized waveform 
undergoes disintegration, transforming into a quasi-linear nonstationary extended state. This state keeps expanding 
until it eventually hits edge absorbers, thus completely disappearing. Naturally, the expansion occurs faster in the 
temporal (alias longitudinal) direction, where it is not impeded by any potential structure. Figure 01 illustrates the 
disintegration of the localized state (the one from Figure ^ , following gradual decrease of e along the propagation 
direction. It is relevant to stress that the disintegration of the soliton at small e is inevitable, as Eq. with e = 
has no 2D soliton solution, unlike the Townes soliton, which would be a solution for the equation with e = and 
reverse sign in front of the u^x term. 

Detecting the delocalization transition of the semi-gap solitons at critical values of the nonlinear coefficient (or 
rescaled norm) and/or the strength of the periodic potential can be used to locate the lower border of their existence 
region in the parameter space (TV, e). Actually, the transition from the localized state to the extended one is quite steep 
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FIG. 3: A stable localized solution of Eq. (0) with e ~ 2.0, produced from the initial Gaussian pulse H13() 
with the initial norm A'o — vr. The amplitude, inverse width along the temporal direction, and norm of this 
localized state are A — 0.74, a = 0.56, and N = 2.29, respectively. Note that 27% of the initial norm was 
lost with emitted radiation in the course of the evolution. Notice that the soliton features no extended tail. 




FIG. 4: Dynamical disintegration of the localized state shown in Fig. as a result of a gradual decrease of the 
strength of the periodic potential along the propagation distance, so as e{z) — £o{l — z/zend), with eg = 2 and 
^cnd — 1800. Cross sections of the wave profile are displayed: (a) along the transverse coordinate x; (b) in 
the temporal direction t. The delocalization occurs around z — 650, at e close to a critical value, Sct — 1.3. 

(see Fig. which allows quite accurate determination of the critical values of the parameters. However, delineating 
the full existence region of the semi-gap solitons is a harder problem, as the limiting effect at large values of the norm, 
which determines the upper border, is splitting of the soliton (see below), rather than collapse (singularity formation) 
in the case of lattice gap solitons with the self- focusing nonlinearity Q . Precise shapes of the existence and stability 
domains of lattice gap solitons can be rather complex, as shown in Ref. for the case of saturable nonlinearity. 

Accurate determination of the full stability borders for the solitons in the present model, going beyond the use of 
the VK criterion and collection of typical examples of direct simulations, will be a subject of separate work. 

Finally, Eq. (|2Jl features obvious Galilean invariance in the longitudinal it) direction, which makes it possible to 
generate a boosted soliton Uc, with an arbitrary inverse- velocity shift c, from a given soliton u, as 

Uc{z, t) = u{z, t - cz)e-<'='/^)^-'"*. 

The use of such two pulses with the ci 7^ C2 makes it possible to study collisions between the moving solitons, which, 
however, should be a subject of a separate work. 
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FIG. 5: The amplitude of the 2D semi-gap sohton abruptly decays at the point of the delocalization 
transition, which was displayed in Fig. ^ The present figure illustrates the steepness of the transi- 
tion at e ~ 1.3, with the slowly decreasing strength of the periodic potential, starting from e — 2.0. 

IV. EFFECTS OF NONLINEARITY MODULATION 

As the nonlinearity is a key factor necessary for the existence of solitons, in this section we address response of 
the 2D solitons to variation of the nonlinearity strength along the propagation coordinate, z. In nonlinear optics, a 
variable nonlinearity coefficient can be created by dint of different physical mechanisms, such as variation of a dopant 
concentration, optically controlled photorefraction, or simply by using a variable thickness of the planar waveguide. 

Thus, we replace Eq. (j^J by its variant which includes a variable nonlinear coefficient x(2;), 

iuz - ^utt + ^Uxx +ecos{2x)u - x{z)\u\'^u = 0. (14) 

Then, we follow the transformation of the soliton with slow increase or decrease of x{z)- FigureEldisplays the evolution 
of the soliton's parameters as the nonlinearity coefficient gradually increases ten-fold. 
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FIG. 6: The evolution of the norm {N), amplitude (A), and inverse temporal width (a) resulting from the slow linear 
increase of the nonhnear coefficient x in Eq. (0), x{z) — 1+jz/zend, with 7 = 9 and < z < ^cnd = 2000. Note signifi- 
cant decrease of the norm N, opposed by little variation of A and a. The initial localized state is the same as in Fig. O 

In this case. Fig. [7| shows that the final waveform, corresponding to x = 10, has not changed notably compared 
to initial one (see Fig. which implies that the increase of the nonlinearity is countered by the loss of the norm. 
The excess norm is shed off with linear waves which are absorbed on the domain boundaries. It appears that the 
shape of the localized wave, being weakly sensitive to the value of the norm (hence, to the strength of the nonlinearity 



8 



too), is actually fixed by the strength e of the periodic potential. This shape- invariance property of the mixed- type 
solitons is very different from what is manifested by both ordinary solitons and gap solitons per se, whose shapes are 
particularly sensitive to the strength of the nonlinearity, at a fixed amplitude of the periodic potential. 




FIG. 7: The final shape of the soliton from Fig. 13 established after by slow in- 

crease of the nonlinearity by a factor of 10, as illustrated in the previous figure. Note 
that an extended tail attached to the soliton does not appear in this case either. 

If the coefficient of nonlinearity slowly decreases, as x(-z) = 1 — z/zond with 2;ond — 2000, the disintegration of the 
localized state is observed when falls to the level of x — 0.5, resembling the picture in Fig. 0] Thus, we conclude 
that the strengths of both the periodic potential and nonlinearity must exceed some critical values in order to sustain 
the localized states. In this respect, the multidimensional semi-gap solitons resemble regular gap solitons in lattice 
potentials [l7llT9l|. 

One of characteristic features of ordinary solitons in ID nonintegrable systems is splitting of the soliton when the 
GVD coefficient 32] or nonlinearity js^ is abruptly changed. In the present model, one can observe a similar effect 
for 2D solitons. Figure |HI displays an example of the splitting of a soliton into two fragments, which then separate 
along the free direction t, after the nonlinearity coefficient was suddenly increased by an order of magnitude. 




t 

FIG. 8: Splitting of the soliton from Fig. [21along the free direction t, following the rapid increase of the nonlinear coeffi- 
cient by a factor of 10, which is performed by setting x(z) ~ 1 -f 7 tanh(4z/zcnd) in Eq. (fnji . with 7 = 9 and Zcnd = 10. 



V. THE THREE-DIMENSIONAL CASE 



The 3D version of Eq. has a straightforward form, 

iuz ~ ^utt + ^{uxx + Uyy) + £[cos(2a:;) + cos{2y)]u - \u\'^u = 0. (15) 

Similarly to the 2D case, stationary solutions to Eq. (|15|l can be numerically found by using a Gaussian pulse as 
the initial condition and propagating in z. Stability of the 3D solitons was verified by simulating the evolution of a 
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soliton with a random perturbation added to it. As well as in the 2D case, the simulations were performed with the 
absorbers placed at borders of the integration domain (under the condition that the length of the domain was much 
larger than a characteristic size of the soliton). As a result, it was concluded that robust 3D solitons exist as generic 
solutions, without any visible tails attached to them. An example of a stationary 3D soliton, which was found to be 
quite robust in stability simulations, is displayed in Fig. |51 In particular, Fig. |5fb) clearly demonstrates that, as 
well as in the 2D case, the soliton has, roughly, equal sizes in the spatial and temporal directions, i.e., the GVD and 
diffraction modified by the periodic potential play equally important roles in supporting the solitons. 




X 




X 

FIG. 9: A three-dimensional stationary soliton solution to Eq. H15|) with e = 2.0 is shown through its two 
cross sections: (a) perpendicular to the free direction, in the t = plane, and (b) parallel to the free di- 
rection, in the y = plane. The soliton was obtained by direct propagation in z of an initial Gaussian 
with the norm Nq ~ 2ti. The norm of the established soliton is iV = 4.16, i.e., a third of the initial 
norm was lost in the course of the adjustment of the initial pulse to the stationary shape of the 3D soliton. 

Besides the fundamental 3D pulses, such as the one displayed in Fig. El the model can also support 3D solitons 
with embedded vorticity. Analogy with known results for gap-soliton vortices in the 2D lattice models suggests 
that the vortex solitons too may easily be stable in the present model. However, detailed investigation of the vortex 
solutions, as well as collecting systematic data about the family of the fundamental 3D solitons, requires numerous 
runs of lengthy simulations of the 3D equation, which is beyond the scope of the present paper. 

VI. CONCLUSION 

In this work, we have proposed a new type of the multidimensional model in nonlinear optics. It combines self- 
defocusing nonlinearity and normal group-velocity dispersion with periodic modulation of the local refractive index 
in the one or two transverse directions (in the 2D and 3D models, respectively). Strictly speaking, multidimensional 
(spatiotemporal) solitons cannot exist in media of this type, as the system's spectrum contains no true bandgap. 
Nevertheless, solitons which seem as completely localized ones are predicted by the variational approximation, and 
found in direct simulations. These solitons are solutions of a mixed type, as in the free (longitudinal, alias temporal) 
direction they are regular solitons, while in the transverse direction(s) they are objects of the gap-soliton type (hence 
the solution as a whole was called a semi-gap soliton). The existence of the solitons requires that both the norm of the 
solution (in other words, the nonlinearity strength x) and the strength e of the spatially periodic transverse potential 
exceed certain minimum values, otherwise the pulses decay into linear waves. Actually, the solitons are much more 
sensitive to s than to x- 

The results reported in this paper call for further work, that should be aimed at accurate identification of borders 
of the solitons' stability regions, especially in the 3D model (which requires running very massive simulations), 
investigation of collisions between solitons, that may move freely in the longitudinal direction, and the study of vortex 
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solitons in the 3D case. 
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